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ABSTRACT 



I consider the problem of the evaluation of the orthopositronium decay rate, including 
second order radiative corrections. I present a brief theoretical discussion of this problem and 
a review of the results available in literature, and, then, I analyze some 0(a 2 ) annihilation- 
type contributions, recently computed. 



1. Introduction 
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Positronium is a bound state of an electron and a positron. It is essentially a pure 
quantum electrodynamical system; in fact the effect of strong and weak interactions can 
be safely neglected at the level of accuracy at which we are interested. 

Its study has always been considered an interesting test of the accuracy of QED 
calculations and especially of the formalism used to describe bound states in quantum 
field theory. 

The lifetime of this atom can vary in a quite large range from 10 -10 to s, 



according to its spin state. In fact it can exist in two different spin states: a singlet S — 
state, called parapositronium, and a triplet S = 1 state, called orthopositronium. 

Parapositronium decays mainly into a couple of photons. The triplet state, instead, 
decays electromagnetically into an odd number of photons greater than one, because 
decays into an even number of photons are forbidden by charge conjugation invariance and 
one photon decay would violate the energy momentum conservation. Hence, the channel 
with three photons in the final state is the dominant decay mode for the orthopositronium 
and one can restrict the analysis to this channel at the 10 ppm level of accuracy. 

We will focus our attention to the analysis of orthopositronium decay width. 
Despite of the fact that different precision measurements and theoretical compu- 
tations of this quantity have been performed in the past years, there is still a certain 
indetermination in the results and one cannot exclude a possible discrepancy between 
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theory and experiment. This discrepancy could eventually indicate also a problem in the 
formalism used to study this electromagnetic bound state system. 

The more precise experimental determinations of orthopositronium decay rate, 
obtained by the Michigan University group 3, are the following: 

X e op S = 7.0514 ± 0.0014 fxs' 1 ; A^ s = 7.0482 ± 0.0016 /is" 1 . (1) 

The first measurement was performed by stopping in a gas the positrons produced 
by the (3 decay of a radioactive source, in such a way to produce positronium, and looking 
at the annihilation 7 rays produced by the O-ps decay. The spectrum obtained was then 
fitted to an exponential function, taking care to begin the fit at a late enough time, so 
that the perturbed singlet component was not present anymore. Finally, the values of A 
obtained at different gas densities were extrapolated linearly to zero density. In the second 
Michigan experiment M, instead, a different technique was used and orthopositronium 
was produced in evacuated MgO-lined cavities. 

A more recent Japanese experiment |J found a quite different value 

\™p = 7.0398 ± 0.0025 ± 0.0015 /is" 1 , (2) 

where the first error is statistical and the second systematic. 

There had been also other measurements of the decay width (we can remember 
in particular the value Aq^ = 7.031 ± 0.007 fis~ l obtained by the Mainz experiment 0j 
using an evacuated aluminum cavity), but they were all characterized by higher values of 
the error. 



The theoretical expression for the orthopositronium decay rate, in the approxima- 
tion of considering only the decay into three photons, can be written in the form f\ 
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where the lowest order decay rate A , first obtained by Ore and Powell 0, is given by 

mc 2 2 (tt 2 - 9) 
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7.211169 ± 0.000004 / us~ 1 . 



(4) 



The coefficient A has been computed by many authors H 0, [§ |9| and its most 
accurate determination A = — 10. 286606±0. 000010 has been found in ||. The constants 
^ and of the logarithmic terms have been found respectively in ]7|, |^, |TU| and JTI |. 
Using these known results, one gets: 
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(5) 



We can see from the last expression that there could be a discrepancy between the 
theoretical and the experimental results. In fact, putting B = in eq. ([5|), we obtain a 

'Here we are adopting the convention to collect an explicit power of ( — ) n in front of the contribution of order 
n to the decay width, to be consistent with what is usually done in literature. We would like, however, to stress 
that this convention, justified by the fact that this is the usual natural scale of the relativistic corrections, can be 
misleading. In fact, in this way, the coefficients of the higher order radiative corrections appear unnaturally big. 
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theoretical value at 0(a), with the inclusion of logarithmic terms up to 0(a 3 In 2 a), which 
differs by the experimental results of the Michigan group Q, [1[ by 6.2 a and 9.4 a. To 
reproduce these results one should have a value of the coefficient B ~ 250, that would be 
unnaturally big. On the other hand, B ~ 40 would be enough to reproduce the result of 
the Japanese experiment ||. Hence, up to now, one cannot discriminate whether the so 
called "orthopositronium decay width puzzle" is a theoretical or an experimental problem. 
It is clear, in any case, that a complete 0(a 2 ) calculation is needed. 

Some second order contributions to the decay rate have been already computed by 
many authors, but some other still need to be evaluated. 

Here we focus our attention on some second order annihilation type radiative cor- 



rections that have recently been computed in [|12 



2. Analysis of the bound states in quantum field theory 

The study of a bound state like positronium is a complex problem that can be 
faced only using difficult techniques of quantum field theory and making some appropriate 
approximations Ul3| . 

The first simplification of the problem can be obtained considering positronium 
as a two body system. The development of a consistent relativistic two body formalism 



for bound state calculations is mainly due to the works of Schwinger and Bethe and 
Salpeter |15| . 



To study the bound states of a two body system, we have, essentially, to look 
for the poles of a four point Green function, that, for fermions, is defined in the usual 
quantum field theory way as Q: 

G(x l ,x 2 ,x 3 ,x A ) = (0|T (V^iM^M^M^)) |0). (6) 
This Green function must obey the Bethe-Salpeter equation 

G(xi,x 2 ,x 3 ,x i ) = (xi,x 3 ) (x 2 ,x 4 ) + 

+ S?M x *(*,,*,,„,„) G (»,,*,, 

where Sp is the Feynman propagator. 

This equation can also be represented graphically, like in fig. 1 . 

In eq. (|7|) we have denoted by K the sum of all the two particle irreducible graphs, 
that is usually called the Bethe-Salpeter kernel. 

One can use the translational invariance properties of the system and express the 
Fourier transformed Green function in terms of only three variables: 

/ d Xl dx 2 dx 3 dx 4 exp [i ( PlXl + P2 x 2 - p 3 x 3 - p 4 x 4 )] x G(x u x 2 , x 3 , z 4 ) = 

= (2tt) 4 5( 4 ) (p 1 +p 2 -p 3 - p 4 ) G(p, q, P) , (8) 



^Here we are using, for simplicity reasons, a two particle formalism, instead of a particle-antiparticle one. In 
the specific case of positronium, one should of course charge conjugate one of the particles of our system. 

^Notice that, to simplify the notation, here and in the following we are denoting by dx the four dimensional 
integral d 4 x 
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where ^ ^ 

P = 2 (Pi ~ Pa) ! 9 = 2 (P 3 ~ P 4 ) ; P = P 1 + P 2 = P 3 + P 4 • ( 9 ) 
The next step is the introduction of the two-particle (two-antiparticle) physical 



states \i, P) ± of mass Mi and four momentum Pi = (Ei, P), Ei = y P 2 + Mf. 

Using these states, we can define the Bethe-Salpeter wave functions in the following 

way: 

x f p (p) = J d{x 1 -x 2 ) e ipixi -' X2) (0 |T(> (1) 0ri) ^ (2) (x 2 )) \i, P)± e iP ^ +X2) 

Xt >p (p) = J d(x 3 - x 4 ) e- iq - {xz ~ Xi) ± (i,P\T (4> {1) (x 3 ) ^ ] {x A )) |0> e - m ^ (x3+X4) (10) 

One can show that the Green function G(p,q,P) can be written as a sum containing 
the Bethe-Salpeter wave functions x ± an d an d this sum has some discrete singularities 



for specific values of P (for which Pq = ±y P 2 + Mf ). 

Using Bethe-Salpeter equation, we get an equation for the residua of the poles of 
the Green function, that gives: 

xtpip) = s F (ip + p) s F (ip - p) J-0p K(p, s, P)xtA s ) ■ (11) 

The values of P giving non trivial solutions of eq. (O) are the energy-momenta of 
the bound states coupled to the two particle-antiparticle states. 

One usually looks for perturbative solutions of this equation, by writing the kernel 
K as the sum of a piece Kq, whose solution is known, and a remainder AK. 

We can also take advantage of the fact that a QED bound state, like positronium, 
can be formed only if the relative momentum of the e + e~ pair is small, since the binding 
force is weak. The expectation value of the relative momentum in the bound state is of 
the order p ~ a m, where m is the mass of the electron. Hence we are allowed to make a 
non relativistic approximation and to distinguish a "large component" from a small one 
in the solution if) of Dirac equation: 

* = [ X I • (12) 
\ X 

\ po+m /v 

The first order Bethe-Salpeter kernel contains in our case the product of a photon 
propagator D^ v times two Dirac gamma matrices 7^ and 7^ and the only term of this 
product connecting the large-large components is the term D o7o7o- Hence, working in 
the Coulomb gauge, we can choose the following expression for the unperturbed kernel of 
the Bethe-Salpeter equation: 

K (p ,p) = —, ^7o><7o- (13) 

(P ~P) 

In this way, we get for the large-large components the well known Schrodinger 
equation with a coulombic potential. The lowest order expression of the ground state 
wavefunction for the orthopositronium can, therefore, be written as 



a ■ e r 




(14) 
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where 



m = ^ { p2 + 1 2y 5 = \/- ; 7 = — • (15) 
This is the expression that we will use to compute the orthopositronium decay 



rate. 



3. The orthopositronium decay width 

At lowest order, the orthopositronium decay width can be obtained simply consid- 
ering the diagram of fig. 2. 

We can easily see that it is given by 

2M J 2lu\ (2vr) 2u 2 (2vr) 2uj 3 (2tt) eii£2ie3 3 £m 3! 

(16) 

where M is the orthopositronium mass, P = (M, 0, 0, 0) ~ (2m, 0, 0, 0) the rest-frame 
orthopositronium energy-momentum vector and uoi = k® = fcj the energies of the final 
state photons. In eq. flTjD we have denoted by Si and e m , respectively, the polarization 
states of the final photons and of positronium and with Ai the matrix element for the 
graph of fig. 2. 

One recovers, in this way, the well known lowest order value for the decay width, 
already found in || (see eq. (|)). 

At the next order of perturbation theory, we must take into account different classes 
of diagrams, like self-energy and vertex corrections and the radiative corrections given by 
the two graphs of fig. 3, usually denoted as annihilation and binding diagrams. 

The last one is particularly important for many reasons. In fact one can easily 
show that the matrix element M.b for this graph can be written in the form 

M B = M)(l- 3 , (17) 

where M.q is the matrix element of the lowest order graph of fig. 2. 

Hence, it seems that this diagram gives contribution to the decay width, not only 
at order a, but also at order zero. This unusual fact can be explained considering that, 
when we add an additional binding photon to the original lowest order diagram of fig. 2, 
we are considering also the possibility that this binding photon is a coulombic one. On the 
other hand, in the determination of the positronium wave function one has already taken 
into account the exchange of any number of coulombic photons between the electron and 
the positron lines. It is, therefore, clear that one must subtract from the contribution 
of the last graph of fig. 3 the part corresponding to the exchange of a Coulomb photon 
and this correspond to subtract the term of order zero from the expression of eq. ( |TTD . In 
this way, one gets the following result for the matrix element of the "subtracted binding 
diagram" *M' B Q 

^Note that here and in the rest of the paper we write explicitly in the formulas the powers of a /it appearing 
in all the amplitudes; on the contrary, we omit them in the text, with the exception of this line. Note also that 
we will not write any power of a/ir for the two unsubtracted amplitudes Mb and Mab, since they contain terms 
of different order in a /it. 
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a A A' O a A A 

-M B = -3 - M . 

71 71 



The binding diagram is particularly important also because it gives a contribution 
to the decay rate, which is bigger than the 90% of the total 0(a) radiative corrections ||. 



The annihilation diagram gives a contribution |T6| that has the same sign of the binding 
contribution, but is, anyway, smaller than 10% of the total 0(a) correction. The sum of 
all the vertex corrections has the same order of magnitude and the opposite sign, and the 
self energy correction is even smaller [^|, 171]. 

As already said, we must consider also the second order radiative corrections, some 
of which have been already computed. 



A first relevant contribution of this order, first evaluated in |T^] and recently up- 
dated in ||, is given by the sum of the squares of all the first order amplitudes. They 
give a contribution equal to 28.860 ± 0.002 to the coefficient B of eq. (|3|). 

A second contribution |HJ to B, equal to 9.0074 ±0.0009, comes from the radiative 



corrections to the light-light scattering block. The inclusion of the vacuum polarization 
corrections to the first order graphs enhances the value of the coefficient B of 0.964960 ± 
0.000004 as proved in 



Considering the decay channel into five photons, one gets |2T| an additional con- 
tribution to B equal to 0.187 ± 0.011. 

The very important second order relativistic corrections have been studied by differ- 
ent authors. Khriplovich and Milstein [2^j have found a big contribution to the coefficient 
B, equal to 46 ± 3, in agreement with the result of Faustov et al.f23]|. Quite a different 
result have been found with a different approach, in the second paper of fill , by Labelle 



et al., that, using the so called "Non Relativistic Quantum Electrodynamics" have got for 
this contribution the value 24.6. 

Generally speaking, we can write the matrix element for the sum of all the diagrams 
contributing up to 0(a 2 ) in the following way: 



M = M + ^{M B + M a + Mx)+ (M' AB + M AR + M 2 ) + 0(a 3 ) , (19) 



where M\ represents the sum of all the first order amplitudes with the exceptions of the 
first order annihilation diagram , denoted by Ma, and the subtracted binding ampli- 
tude, M B . The second order annihilation type corrections are given by the subtracted 
binding diagram, M AB (fig. 4(A)) and the radiative corrections to the light-light scat- 
tering block, Mar (an example of which is given in fig. 4(B)); M 2 denotes the remaining 
(non-annihilation type) second order amplitudes. 

We have seen that the contribution coming from the first order binding diagram, 
represents more than 90% of the first order radiative corrections; hence it seems reasonable 
to look at second order corrections obtained by graphs containing some additional binding 
photon. 



In a recent paper [T2| , we have examined the contribution to the decay width com- 
ing from the interference between the subtracted annihilation binding graph represented 
in fig. 4(A) and the zero order diagram. In the same paper we have also considered the 
square of the first order annihilation amplitude Ma, contributing at 0(a 2 ) to the decay 
rate, and we have verified the existence of a logarithmically enhanced contribution arising 
from the radiative correction to the light-light scattering block depicted in fig. 4(B). 
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Let's recall briefly how one can evaluate these different contributions (for more 
details look at |I2|). 



The matrix element of the annihilation binding diagram can be written as (see fig. 

4(A)): 

M { T B X) = --^-T^G^ p . (20) 

In the previous formula the tetravector G^ e ' p describes the transition of the heavy 
photon to the three real ones and we have denoted with e = (si, e 2 , e 3 ) the set of the three 
polarizations of these photons (e$ = ±1). The symbol represents the 0(a) correction 
to the annihilation current 4- vector of the positronium in the polarization state e m . 

One can easily see that contains a double integral, over the variables p and 
k, a trace of 7 matrices, including also the orthopositronium wave function ^/^ m \p), and 
one photon propagator that we can write as — % 

The A MI/ tensor depends on the gauge we use. The choice of the gauge is subtle 
when dealing with bound state problems. It has been discussed elsewhere (see for instance 
the last paper of [7J) that the Coulomb gauge is the most natural for calculations in 
positronium. However, covariant gauges are simpler for computing radiative corrections, 
and, among them, the Fried- Yennie (FY) gauge 03 is the most convenient, due to its good 
infrared behaviour. We have computed both in the FY gauge and in the Coulomb 
gauge. As expected, the result is the same in both cases, and no gauge correction term 
must be added when using the FY gauge. 

Let's report the basic steps of the calculation in the FY gauge, for the analogous 
computation in the Coulomb gauge we refer the interested reader to |12 . 

In the FY gauge we have A^ v = g^ v + 2-^. To perform the computation in this 
gauge, we have splitted the trace entering T^" 1 -* into two pieces, one remaining non-singular 
at k = and another one containing the contribution of the coulombic photon. Formally 
we have used the following equality: 

Tr ^(k) = Tr ^(0) + (Tr „ vp (k) - Tr ^(0)) . (21) 

The first term gives a contribution to the matrix element proportional to the 0(a) 
annihilation amplitude ^M-a 

M AB>1 = -M A (l - 3-) . (22) 



TT \ TT 



The second term of eq. ( pT|) is infrared finite and we can safely put p = in the 
loop integral, introducing an error of order 0(a 2 ). The ultraviolet divergence introduced 
by this second term can be regulated either with dimensional regularization or with the 
use of a cut-off. 

One must also add the contribution of the "annihilation vertex" counterterm, that 
cancels the divergence coming from the term (Tr fiup (k) — Tr Mi , p (0)) as explicitly proved 
in ||r^| . The sum of the contributions of the second term in eq. (^) and of the counterterm 
gives: M A b,2 = +~M A . 

Hence we have found the following expression for the unsubtracted annihilation 
binding diagram in the FY gauge: 

M AB = (l-2-)-M A . (23) 

V IT J TT 
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We must subtract the lowest order contribution - Ai A corresponding to the ex- 
change of a coulombic binding photon (as explained in the case of the 0(a) binding 
correction) and consider the interference of this subtracted annihilation binding diagram 
with the zero order graph. In this way we have got the following 0(a 2 ) contribution to 
the orthopositronium decay width: 



ci' .. „ „ /a N 2 



r AB = -2-r A = 1.6281 - r , (24) 

7T \7T/ 

where in the last equation we have used the numerically improved value of the lowest order 
annihilation width, that can be found in the first paper of [p0[| : Ta = —0.81406 (oc/tt) Tq . 



This result is in very good agreement with the estimate, based on factorization arguments, 



of this correction, that was made by Karshenboim in [25 



In |T2[ we have computed also the contribution to the decay width coming from the 



square of the order 0(a) annihilation amplitude (see fig. 3(A)). This square contains the 
integral over the phase space of the final photons of the product of two tensors G £1,e2,es . 
To evaluate this quantity we have used the following relations (that can be found 



in [|6| and in |7j) 



U\ UJ2 ^3 



and 



(I 'M'Vv >^3 5 (4) (p - kl -k 2 -k 3 )=8 vr 2 m 2 J dv x dv 2 du 3 5(2 -v x -v 2 - v 3 ) (25) 
- GfG {£)p = 2 6 a 4 [#(234) + #(324) + #(423)] , (26) 



where v i = ^ = L-U and R are complicated functions of Ui computed in |27j and reported 
in 0. 

Using eq. (|25|) and eq. (|2~6|) we reduced the evaluation of the square of the anni- 
hilation amplitude to the calculation of a numerical integral. We have got the following 
numerical result: 

T A 2 = (0.17021 ± 0.00010) T . (27) 



Let's notice that this result, like the annihilation-binding contribution of eq. ( 24]) 
(and differently from the lowest order annihilation contribution of 0(a)) has the right 
sign to reduce the eventual discrepancy between theory and experiment. Nevertheless, 
the absolute value of these corrections to the decay width is quite small and they are 
manifestly far from solving this discrepancy. If the "orthopositronium problem" has to 
be solved by this kind of perturbation theory, larger contributions must be searched in 
other classes of diagrams. 



Finally, let' s remember that in [T2| we have also considered the radiative correc- 
tion to the light-light scattering block given by the graph of fig. 4(B). It generates a 
logarithmically enhanced term that gives a contribution to the decay width proportional 
to a 2 ln(a)r . 

To recover this result one can consider that there are two regions of the loop 
momenta space giving the main contribution: one corresponding to /, k ~ m and another 
one with fermion momenta almost on mass-shell, that is with /, k ~ am. The logarithmic 
term is produced in this second region. In the analysis of this region one can put almost 
everywhere in the integral k and 1 to zero. The integral over / would be ultraviolet 
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divergent, but, according to the considerations we have just made, we can introduce the 
ultraviolet cutoff m, and get the correct infrared logarithmic term. We have found a 
contribution to the decay width that can be written as 



Ar 



AR 



-a 2 In — r n 



a 



(2f 



in agreement with the results of [19], where the all set of these radiative corrections has 
been computed, and with the ones of |TT1 and of the first paper of [171. 
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G 




Figure 1: Graphical representation of the Bethe-Salpeter equation. See the text for the 
meaning of the symbols. 
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Figure 2: Graph determining the orthopositronium decay width at lowest order. 
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Figure 3: Two different corrections of 0(a):the annihilation graph and the binding dia- 
gram. 
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Figure 4: Two different kinds of corrections to the annihilation graph. (A) is the vertex 
correction, (B) represents the insertion of a photon into the light-light scattering block, 
which generates a logarithmically enhanced contribution. The direction of the fermions 
in the final loops is clockwise for all the two diagrams. 
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